AD-A221  317 


WRDC-TM-90  -  300 

MERC2PLOT3D  User's  Manual 
W.  Z.  STRANG 


AERODYNAMIC  METHODS  GROUP 
AERODYNAMICS  &  AIRFRAME  BRANCH 
AEROMECHANICS  DIVISION 


JANUARY  1990 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION 
IS  UNLIMITED. 


FLIGHT  DYNAMICS  LABORATORY 

WRIGT  RESEARCH  AND 
DEVELOPMENT  CENTER 

WR IGHT-P ATTERSON  AFB ,  OH  45433 


90  05 


09  129 


FORWARD 


This  effort  was  accomplished  under  Work  Unit  240410A1,  entitled 
"Aerodynamic  Design  and  Analysis  Methods."  The  effort  covers  work 
performed  during  January  1990. 
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The  Euler  and  Navier-Stokes  flow  solvers  of  computational  fluid 
dynamics' /fCFD)  generate  so  much  output  data  that  one  needs  interactive 
computer  graphics  to  fully  and  quickly  digest  the  results.  Multi-block 
solutions  can  be  even  more  difficult  to  interpret  due  to  the  patching 
required  at  the  block  abutments.  Fortunately,  NASA  Ames  developed  an 
excellent  graphics  post-processor  of  CFD  solutions,  PL0T3D  meeting 

the  requirements  of  multi-block  solutions.  Thus,  all  that  is  required 
to  inspect  multi-block  flow  solutions  (when  one  has  PL0T3D  and  a 
suitable  workstation)  is  to  convert  the  output  from  the  CFD  solver  to 
natch  the  input  requirements  of  PL0T3D. 

MERC2 PL0T3 D  converts  MERCURY restart  files  into  the  proper 
input  data  for  PLOT3D.  —  '  i  i  _  r  •'  .  I 

-■''r  '  '  ‘  1  ,  J 

Running  MERC2PLOT3D  v 

MERC2PLOT3D.  FOR  is  run  in  batch  mode  on  a  Cray  computer.  The  user  'A. 
must  modify  the  JCL  of  MERC2PLOT3D. JOB  to  access  the  MERCURY  restart  ' 

file  on  the  Cray  and  fetch  the  appropriate  connectivity  file  [2]  from 
the  Vax  front-end.  When  changing  file  names,  change  only  the  PDN ,  ID 
and  ED  parameters;  leave  the  DN  parameter  alone.  PL0T3D  requires  the 
flow  solution  data  and  the  grid  data  to  reside  in  different  datasets. 

The  user  must  also  specify  the  names  of  these  two  datasets  which  will  be 
saved  on  the  Cray  upon  completion  of  the  MERC2PLOT3D  run.  See  Fig.  1. 

Following  the  JCL  in  MERC2PLOT3D. JOB  is  FORTRAN  unit  5  (or  $IN) 
which  MERC2PLOT3D  reads  as  input.  See  Fig.  2.  The  required  input 
should  be  self-explanatory  except  for  the  question  regarding  the  version 
of  MERCURY  that  created  the  restart  file.  On  27  December  1989,  the 
MERCURY  restart  file  structure  was  changed  in  order  to  reduce  the 
restart  file  size  by  11%.  Thus,  edit  MERCURY. FOR  to  determine  the  date 
of  the  last  revision.  This  date  occurs  on  line  eight.  If  that  date  is 
27  December  1989,  or  later,  the  correct  response  is:  0.  Otherwise,  your 
restart  file  was  created  by  an  old  version  of  MERCURY  and  the  correct 
response  is:  1. 


JOB , JN-PIX , MFL , T-500 , SSD-IOOOO . 

ACCOUNT , AC-D840262 , APW-SIRBAUGH . 

OPTION, STAT-OFF . 

* . 

*.************************************************** 
*.  ACCESS  THE  MERCURY  RESTART  FILE  ON  THE  CRAY 
*.********** **************************************** 
* . 

ACCESS , DN-REST , PDN-CYLNDRAOR , ID=GREG . 

* . 

*. ************************************************** 
* .  FETCH  THE  CONNECTIVITY  FILE  FROM  THE  FRONT-END 
*.  VAX  ACCOUNT 

*.************************************************** 
* . 

FETCH , DN-CAPRI , TEXT- ’ CYLNDRCN . DAT 9 . 

* . 

*.************************************************** 
*.  LEAVE  THE  FOLLOWING  EIGHT  JCL  LINES  ALONE 
*  # ************ ************************************** 
* . 

FETCH , DN=COMET , TEXT= ’ MERC2PLOT3D . FOR ’ . 

ASSIGN , DN=REST , A=FT01 , LM=200000 . 

ASSIGN , DN=QOUT , A=FT02 , LM=200000 . 

ASSIGN , DN=XOUT , A=FT03 , LM=200000 . 

UPDATE , P=0 , I=CAPRI : COMET . 

CFT77,I=$CPL, DEBUG. 

SEGLDR , L=0 , CMD= ’ ABS=PIX ’ . 

PIX. 


*  .  ************************************************* 
*.  SAVE  THE  Q  AND  X  PL0T3D  FILES  ON  THE  CRAY 
* . ************************************************* 
* . 

SAVE , DN=QOUT , PDN=CYLNDR_AO , ID=QPL0T3D . 

SAVE , DN=XOUT , PDN=CYLNDR_AO , ID=XPL0T3D . 

* . 

AUDIT. 

EXIT. 

DUMP JOB. 

DEBUG. 

/EOF 


Figure  1 


I/O  DEVICE  (0=SSD,  1=DISKS) 

0 

RESTART  FILE  OF  A  MERCURY  VERSION  PRIOR  TO  27  DEC  89  (0=NO,  1=YES) 
1 

HOW  MANY  BLOCKS  ARE  TO  BE  OUTPUT? 

4 

OUTPUT  WHAT  BLOCKS? 

4 

1 

3 

2 

MACH  NO.  ALPHA  BETA  DT 
0.20  0.0  0.0  0. 

RATIO  OF  SPECIFIC  HEATS 
1.4 


Figure  2 


MERC2PLOT3D  allows  one  to  output  any  number  of  blocks  (you  don't 
have  to  output  all  the  blocks)  in  any  order.  For  example,  consider  a 
multi-block  solution  and  grid  of  four  blocks.  The  user  can  specify 
MERC2 PL0T3 D  to  order  the  blocks  on  output  such  that  block  four  is  the 
first  block  output,  block  one  is  the  second  to  be  output,  block  three  is 
the  third  and  block  two  is  last  by  listing  the  blocks  to  be  output  as: 

4 

1 

3 

2 

See  Fig.  2.  Reordering  the  blocks  upon  output  is  helpful  when  dealing 
with  complex  configurations  where  the  solid  surfaces  defining  the 
vehicle  of  interest  are  scattered  throughout  the  multi-block  grid  in  a, 
usually,  random  fashion. 

MERC2PLOT3D  creates  an  IBLANK  [1]  array  for  each  grid  block. 

PL0T3D  uses  the  IBLANK  values  to  properly  patch  flow  features  (such  as 
streamlines)  across  block  boundaries.  Thus,  to  properly  read  the 
datasets  created  by  MERC2PLOT3D  into  PLOT3D,  use  the  "read/mgrid/blank" 

command . 

Lastly,  MERC2PLOT3D  uses  the  Cray  I/O  devices.  Further,  it  expects 
those  devices  to  have  certain  device  names.  Therefore,  if  you  are  not 
using  the  ASD  Cray  you  must  edit  MERC2PLOT3D. FOR  and  search  for 
.SUBROUTINE  SCRATCH.  In  this  subroutine  you  will  find  the  ASD  Cray 
device  names  spelled  out.  Simply  swap  those  names  with  the  device  names 
at  your  installation.  If  you  are  using  the  ASD  Cray,  be  aware  that  if 
the  SSD  is  down  and  you  are  forced  to  use  the  disks,  and  one  or  more  of 
the  disks  are  down,  MERC2PLOT3D  will  bomb.  If  this  does  happen,  inspect 
the  dayfile.  In  it  you  will  see  the  device  name  of  the  disk  that  is 
down.  Then  edit  MERC2 PL0T3 D ,  find  SUBROUTINE  SCRATCH,  and  swap  the  name 
of  the  down  disk  for  one  that  is  operational. 
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